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ABSTRACT 

In  this  report  a  three  dimensional  finite  element  (F.E.)  model  has  been  developed  for  a 
structural  detail  in  an  F-111  lower  wing  skin.  The  location  of  interest  is  a  fuel  -  flow 
groove  in  the  lower  wing  skin  where  cracking  had  occurred  in  service  on  aircraft  A8- 
145.  Detailed  models  were  developed  for  (i)  xm-cracked  structure;  (ii)  the  cracked 
structure  and  (iii)  the  repaired  structure,  using  a  bonded  composite  patch  for  the 
repair.  The  objective  of  this  work  is  to  validate  the  design  analysis  used  by  tiie  RAAF, 
using  an  independent  approach  for  the  stress  analysis.  The  F.E.  model  has  been 
validated  using  strain  gauge  restdts  from  a  full  scale  test  wing.  The  results  of  the  F.E. 
analysis  are  shown  to  compare  favourably  with  closed  form  solutions  used  by  the 
RAAF  (  RAAF  Engineering  Standard  C5033  )  in  the  original  design  of  the  repair.  Thus 
the  present  work  provides  a  basis  for  confidence  in  the  design  procedures  contained  in 
RAAF  Engineering  Standard  C5033. 
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Finite  Element  Analysis  of  an  F-111  Lower 
Wing  Skin  Fatigue  Crack  Repair 


Executive  Summary 


Following  the  discovery  of  a  48  mm  fatigue  crack  in  an  Australian  F-111  aircraft,  the 
RAAF  implemented  a  boron/ epoxy  composite  repair,  in  which  the  design  was  based 
on  closed  form  solutions  using  RAAF  Engineering  Standard  C5033.  Since  primary 
structure  was  involved,  the  RAAF  requested  that  the  design  analysis  be  validated  by 
an  independent  method  and  a  comprehensive  structural  testing  program  involving 
both  static  and  fatigue  loading.  The  objective  of  the  present  work  being  the  validation 
of  the  design  analysis  using  fkiite  element  analysis  (F.E.)  as  the  independent  approach 
for  the  stress  analysis. 

In  this  report  a  three  dimensional  F.E.  model  has  been  developed  for  a  structural  detail 
in  the  F-111  lower  wing  skin  which  contains  a  fuel  flow  groove  in  which  cracking  has 
occurred.  A  detailed  F.E.  model  has  been  developed  for  the  un-cracked  structure  and 
validated  by  a  static  strain  survey  on  a  full  scale  test  wing.  This  model  has  been 
extended  to  include  a  48  mm  crack,  which  provides  a  basis  for  residual  strength 
calctilations  to  be  made.  Finally,  the  adhesive  and  boron/ epoxy  composite  repair  has 
been  added  to  the  model,  and  has  also  been  validated  by  strain  gauge  data. 

The  results  obtained  from  these  models  compare  favourably  with  closed  form 
solutions  used  in  the  original  RAAF  design,  (using  C5033)  and  the  present  work 
provides  a  basis  for  the  confidence  in  the  design  analysis  of  C5033.  Furthermore,  the 
results  given  here  show  that  the  static  strength  is  restored  and  the  crack  growth  rate  is 
substantially  reduced.  The  results  from  this  analysis,  together  with  the  results  from  the 
extensive  structural  testing  program,  have  substantiated  this  repair  for  the  expected 
service  life  of  the  F-111. 
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Notation 


A  Coefficient  in  equation  (10) 

B  Coefficient  in  equation  (10) 

D  Bending  stiffness 

E  Young's  modulus 

El  Yovmg's  modulus  for  skin 

Eo  Young's  modulus  for  patch 

Ec  Transverse  tensile  modulus  for  adhesive 

E'c  Effective  transverse  tensile  modulus  of  the  adhesive 

Ki  Mode  1  stress  intensity  factor 

Kiup  Mode  1  stress  intensity  factor  unpatched 
Kip  Mode  1  stress  intensity  factor  patched 

Kic  Fracture  toughness 

Ga  Adhesive  shear  modulus 

L  Constant 

N  Cycles 

a  Half  crack  length 

b  Half  panel  width 

1  Crack  tip  element  size 

ta  Thickness  of  adhesive 

ti  Thickness  of  skin 

to  Thickness  of  Boron  patch 

V  Crack  opening  displacement 

Wo  Transverse  displacement 

X  Transfer  length 

a  Direct  stress 

CTg  Direct  stress  in  Boron  patch 

V  Poisson's  ratio 

Y  E  Elastic  limit  for  adhesive  shear  strain 

Y  „  Ultimate  failure  shear  strain  for  adhesive 

r  Shear  stress 

Elastic  limit  for  adhesive  shear  stress 

Equivalent  elastic  failure  stress 
r  Ultimate  failxire  for  adhesive  shear  stress 

P 

Adhesive  shear  stress 
Load  attraction  factor 
Inclusion  factor  for  thermal  stresses 
X  Coefficient  in  equation(l  0) 
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1.  Introduction 

A  finite  element  (F.E.)  analysis  has  been  carried  out  on  a  RAAF  designed  boron/ epoxy 
repair  to  an  F-111  in  which  cracking  has  occurred.  Specifically,  the  crack  was  48  mm, 
tip  to  tip,  nmning  in  a  chordwise  direction  at  Forward  Auxiliary  Spar  Station  (FASS) 
281.28  in  the  lower  skin  directly  beneath  the  forward  auxiliary  spar. 

The  original  design  of  the  repair,  ref  (1),  was  according  to  the  RAAF  Engineering 
Standard  C5033,  ref  (2)  and  is  based  on  equations  with  closed  form  solutions.  The 
objective  of  the  current  work  is  to  validate  the  design  using  an  independent  approach, 
in  this  case  the  Finite  Element  (F.E.)  method.  The  repair  is  also  being  substantiated 
through  an  extensive  experimental  program,  however  this  report  focuses  on  the  F.E. 
aspects. 

The  calibration  of  the  baseline  F.E.  model  used  in  this  report  has  been  described  in  ref 
(3).  Work  presented  here  involves  the  crack  modelling,  evaluation  and  validation  of 
stress  intensity  factors,  and  extension  of  the  analysis  to  model  the  patch  and  adhesive. 
Finally,  a  comparison  will  be  made  between  the  RAAF  Engineering  Standard  and  F.E. 
results. 

Consider  now  the  geometry  of  the  wing.  The  shaded  region  in  Figure  1  (a)  shows  the 
area  of  the  lower  port  wing  skin  detailed  in  Figtue  1  (b).  Shown  in  Figure  1  (b)  are  the 
fuel  flow  holes  in  the  forward  auxiliary  spar  together  with  the  crack  and  patch  outline 
on  the  outboard  location. 

The  crack  in  the  wing  skin  runs  in  a  chordwise  direction  and  ocaurs  in  a  region  in 
which  material  has  been  removed  from  the  integral  stiffener  in  the  skin  to  allow  fuel  to 
move  freely  in  order  to  maximise  fuel  usage.  In  the  design  of  the  wing,  side  stiffeners 
have  been  used  in  an  attempt  to  restore  structural  strength. 
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Inspection  areas 


(b)  Interior  view  of  shaded  area  indicated  in  (a) 


Figure  1  Geometry  ofF-111  lower  wing  skin. 

A  section  view  of  this  area  is  shown  in  Figure  2,  in  which  the  y  axis  is  in  tiie  chordwise 
direction,  the  x  axis  is  in  the  spanwise  direction  and  z  axis  is  in  the  thickness  direction. 
The  crack  is  shaded  and  shown  in  the  yz  plane. 
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Figure  2  Detail  section  view  of  crack  and  fuel  flow  hole  area. 


2.  Modelling 


2.1  Model  description 

Firstly,  it  is  necessary  to  consider  the  geometry  of  the  region  in  which  the  fatigue  crack 
has  been  foimd.  The  crack  runs  in  a  chordwise  direction  and  occurs  in  a  region  in 
which  material  has  been  removed  from  the  integral  stiffener  in  the  skin,  as  shown  in 
Figure  2.  This  was  to  allow  fuel  to  move  freely  and  hence  maximise  fuel  usage.  Side 
stiffeners  which  can  be  seen  in  Figure  lb  and  Figure  2,  were  included  in  the  structure 
in  an  attempt  to  restore  the  structural  integrity. 

Since  the  geometry  in  this  area  is  symmetric  about  the  chordline  which  contains  the 
crack,  it  has  been  possible  to  use  an  F.E.  'half  model'.  This  is  particularly  convenient 
since  it  reduces  the  size  of  the  problem  considerably.  Given  that  the  size  of  the  patch  is 
468  mm  by  324  mm,  die  baseline  model  extends  in  the  chordwise  direction  from  the 
front  spar  to  the  centre  spar.  In  the  spanwise  direction  the  model  extends  beyond  the 
end  of  the  patch  to  a  transverse  member,  which  runs  approximately  chordwise.  The 
dimensions  of  the  baseline  model  are  shown  in  Figure  3. 
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_ >11^ _ 2  mm 


Figure  3  Overall  dimensions  of  baseline  F.E.  model 

As  shown  in  Figure  2,  the  bottom  half  of  the  spar  is  included  in  the  model,  although  in 
this  view  it  has  been  translated  vertically  to  allow  the  geometry  to  be  seen  in  the 
region  of  the  crack.  Note  that  4  radii  are  necessary  to  define  the  geometry  together 
with  the  skin  thickness. 

For  the  repair  substantiation  it  has  been  necessary  to  base  the  F.E.  models  on  the 
engineering  drawings  and  on  minimum  limits  given  the  tolerances.  The  minimum  skin 
thickness  used  is  3.6  mm  and  the  relevant  drawings  are  given  in  ref  (4).  It  is  interesting 
to  note  that  die  skin  thickness  on  a  USAF  F-lllG  cracked  at  the  same  location  was  3.3 
mm.  The  inclusion  of  the  spar  in  the  analysis  is  important  since  the  spar  carries  a 
considerable  part  of  the  load,  and  the  spar  is  connected  to  the  skin  by  interference  fit 
'Taper  lok'  bolts.  Modelling  of  the  spar  also  involves  the  spar  being  separated  by  a 
very  small  vertical  gap  from  the  skin.  This  allows  connections  to  be  made  between 
nodes  on  the  spar  and  skin,  corresponding  to  bolt  locations.  The  importance  of  the 
inclusion  of  bolts  in  the  analysis  can  be  seen  in  deformation  pattern  shown  in  Figure  4. 
Bolt  locations  have  been  fotmd  to  influence  the  secondary  bending  in  this  area.  This  is 
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also  the  location  in  which  fatigue  cracking  occurred.  To  prevent  'merging'  of  the  skin 
into  the  spar  it  has  been  necessary  to  use  gap  elements. 


Figure  4  Side  elevation  of  skin  &  spar  flange.  Deformation  of  skin  between  bolt  centres  is  also 
illustrated. 

Table  1  Material  properties  of  the  spar  and  skin,  ref(5).(  specification  QQ-A-250/4(r851)) 


Material 

E  (MPa) 

G  (MPa) 

V 

Yield  strength 
(MPa) 

Ultimate  tensile 
strength  (MPa) 

2024-T851 

72397. 

27217. 

.33 

407. 

448. 

2.2  Assumptions 

The  following  assmnptions  are  made  for  the  geometry  and  loading  of  the  model. 

1.  Loading  is  in  the  spanwise  direction  only,  and  is  taken  to  be  a  uniform  pressure. 

2.  The  skin  thickness  in  this  model  is  taken  to  be  uniform,  except  for  the  integral 
stiffeners  and  spar  landing. 

3.  The  skin  material  is  taken  to  be  isotropic,  and  whose  properties  are  shown  in 
Table  1. 

4.  Simple  support  conditioirs  have  been  assumed  along  the  edge  of  the  wing  skin 
which  joins  the  centre  and  forward  spars. 

5.  Simple  supports  are  applied  to  the  top  of  the  auxiliary  spar  web. 

6.  Simple  support  conditions  are  assximed  at  die  loaded  end,  which  corresponds  to 
a  transverse  member. 

7.  No  support  conditions  are  applied  to  the  model  in  the  plane  of  the  crack. 
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2.3  Geometry 

The  overall  dimensions  of  the  F.E.  model  are  shown  in  Figure  3.  This  shows  both  the 
plan  view  of  the  half  model,  and  the  dimensions  of  the  spar.  The  definition  of  the  radii 
Rl,  R2,  R3,  and  R4  are  shown  in  Figure  2.  The  drawings  specify  all  diese  radii  to  be  9.6 
mm. 


2.4  Elements 


In  this  analysis  the  use  of  the  prism  element  which  can  result  in  over  stiffening  of  the 
model,  has  been  minimised.  Specifically  for  the  PAFEC  program  the  elements  used 
are: 

1.  20  noded  bricks 

a)  37110  isotropic,  used  to  model  the  skin  and  spar 

b)  37115  orthotropic,  used  to  model  the  adhesive  and  boron 

2.  16  noded  prisms 

a)  37210  isotropic,  used  for  the  skin 

b)  37215  orthotropic,  used  to  model  the  adhesive  and  boron  patch 

3.  Gap  elements 


2.5  Restraints 

Figure  5  shows  the  restraints  applied  to  the  model.  The  directions  X,  Y,  Z  have  already 
been  defined  in  Figure  2.  Specifically  the  following  restraints  are  applied: 

1.  Skin 

a)  Simply  supported  in  the  Z  direction  applied  along  the  sides 
corresponding  to  the  front  and  centre  spars. 

b)  Simply  supported  restraints  in  the  Z  direction,  applied  at  the  loaded  end 
of  the  skin  corresponding  to  a  transverse  member. 

c)  Restraint  in  the  X  direction  corresponding  to  the  Y-Z  plane  for  which 
X=0. 

d)  There  is  one  restraint  in  the  Y  direction  as  shown  in  Figure  5. 
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2.  Spar 

a)  Simply  supported  restraints  in  the  Z  direction  applied  along  the  top  of 
the  spar  web  as  shown  in  Figure  5. 

b)  Restraint  in  X  direction  corresponding  to  the  Y-Z  plane  for  which  X=0. 


Figure  5.  Restraints  applied  to  spar/skin  model 

2.6  Loading 

2.6.1  Flight  loads 

As  stated  in  section  2.1,  it  was  decided  to  base  the  repair  substantiation  on  the 
minimum  allowable  skin  thickness  of  3.6  mm  (this  is  a  conservative  approach).  Ref  (3) 
gives  the  spanwise  applied  loading  (193  MPa)  required  for  the  baseline  F.E.  model  to 
correctly  predict  the  outside  surface  strains  for  a  4.2  mm  skin  thickness  wing  imder  the 
Cold  Proof  Load  Test  (CPLT)  loading.  The  CPLT  loading  results  in  a  local  bending 
moment  of  1.60x10^  inch  pounds,  whereas  the  Design  Limit  Load  (DLL)  condition  at 
PASS  281  is  1.82x10^  inch  poimds  (see  ref  (3)).  The  stress  in  this  location  is  known  to  be 
a  direct  function  of  the  applied  wing  bending  moment.  The  flight  load  to  correspond 
to  a  minimum  skin  thickness  (3.6  mm)  wing  under  DLL  was  therefore  calculated  as 
follows: 

C7  =  193.  C4.2  /3.6;.  (1.82  /1.60;  =  256MPa 

The  bolt  locations  shown  in  Figure  4  were  fotmd  to  have  a  strong  influence  on  the 
results  and  were  adjusted  to  a  distance  slightly  closer  to  the  critical  region  than  the 
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actual  location  to  achieve  the  best  agreement  with  strain  gauge  results.  A  simple 
explanation  is  that  the  clamping  effect  of  the  bolt  extends  a  distance  outside  the  bolt 
diameter. 

2.6.2  Thermal  loads 

The  internal  stresses  and  strains  in  the  repair  system  are  affected  by  the  operating 
temperature  in  two  separate  ways  as  follows: 

(a)  The  mechanical  properties  of  the  cured  adhesive  vary  dramatically  with  operating 
temperature  (see  Table  8).  At  low  temperatures  the  adhesive  is  much  stiffer  and  the 
load  transfer  is  therefore  significantly  different  when  compared  to  room  temperature. 
The  crack  patching  efficiency  is  also  significantly  affected. 

(b)  The  adhesive  bonding  the  patch  to  the  wing  skin  is  cured  at  elevated  temperature 
(82°C).  The  coefficient  of  thermal  expansion  of  fhe  boron  patch  material  is  significantly 
lower  than  that  for  the  aluminium  wing  skin  (4.1x1 0V°C  versus  23.xlO-6/°C).  This 
thermal  mismatch  therefore  results  in  residual  stresses  in  the  system  when  the 
operating  temperature  differs  from  the  bonding  temperature.  The  operating 
temperature  is  usually  lower  than  the  cure  temperature,  so  residual  compression  in  the 
boron,  tension  in  the  aluminium  and  shear  in  the  adhesive  is  expected. 


3.  Stress  and  strain  distribution 

3.1  Strain  distribution 

The  calibration  of  the  baseline  F.E.  model  has  been  carried  out  in  ref  (3).  Briefly  Figure 
6  shows  the  results  of  this  calibration,  where  the  strains  correspond  to  the  location 
on  the  outer  surface  of  the  xmcracked  wing  in  the  critical  region.  Good  agreement  has 
been  obtained  between  experimental  data  and  F.E.  results.  Secondary  bending  is 
shown  by  the  dip  in  the  strain  plot,  Fig.6. 
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Figure  6  Chordwise  variation  of  strain  comparing  FE  results  with  strain  survey. 


The  airworthiness  aspects  require  that  the  repair  is  effective  at  Design  Ultimate  Load 
(DUL)  with  wing  skins  at  the  minimum  thickness  tolerances.  In  this  case  the  minimum 
skin  thickness  is  3.6  mm  and  is  less  than  the  calibration  model  which  has  a  thickness  of 
4.20  mm.  It  has  been  shown  that  the  calibration  F.E.  model  has  given  good  agreement 
with  experimental  results,  and  from  this  a  F.E.  model  has  been  developed  according  to 
the  General  Dynamics  (G.D.)  drawings.  Since  the  only  difference  between  these 
models  is  the  skin  thickness,  it  was  considered  that  this  model  is  also  valid. 


Figure  7  Chordwise  variation  of  strain,  comparing  F.E.  results  for  3.6  mm  skin  with 
calibration  F.E.  model,  at  DLL. 

The  thinner  wing  skin  results  in  higher  stress  levels  and  significantly  greater 
secondary  bending.  The  chordwise  8^  strain  variation  is  shown  in  Figure  7.  Applying 
the  local  bending  moment  and  the  skin  thickness  correction  factors  gives  the  remote 
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stress  at  Design  Limit  Load  (DLL)  as  a*  =  256  MPa  and  applies  to  all  figures  from  this 
point  onwards. 


3.2  Stress  distribution 

From  our  validated  model  the  stress  distribution  can  be  predicted  in  the  area  of  skin  in 
which  cracking  occurs.  The  results  in  Figure  8  show  the  chordwise  variation  of 
spanwise  stress  cr;^ .  (The  location  y  =  0  corresponds  to  the  centreline  of  the  auxiliary 
spar).  On  the  outside  of  the  skin  the  stresses  are  reduced  in  the  vicinity  of  the  spar 
centreline,  but  on  the  inside  very  high  stresses  exist  and  exceed  the  yield  stress  in  a 
localised  area.  (The  yield  stress  for  the  2024-T851  is  400  MPa).  This  explains  why  cracks 
develop  in  this  area,  starting  from  the  inside  and  growing  through  to  the  outside. 
Away  from  this  area  the  stress  field  is  equal  to  the  applied  stress  of  256  MPa. 


Figure  8  Chordwise  variation  of  stresses,  upper  &  lower  surface  of  skin,  at  DLL. 


4.  Deformation  pattern 

A  significant  deformation  pattern  occurs  in  the  area  of  cracking.  As  mentioned 
previously,  bolts  have  been  introduced  in  to  the  solution  in  the  area  of  concern  by 
connecting  specific  nodes  on  the  skin  to  tiie  spar.  The  dotted  line  in  Figure  4  shows  the 
deformation  of  the  skin  as  a  result  of  secondary  bending.  Also  it  indicates  that  some 
bending  of  the  lower  spar  flange  will  occur.  It  has  been  foimd  that  the  bolt  location 
will  have  an  influence  on  the  secondary  bending  since  it  resists  the  shear  between  the 
spar  flange  and  skin.  In  so  doing,  it  resists  the  shear  deformation  and  hence  bending  in 
that  area.  In  fact,  the  calibration  of  the  model  in  ref  (3)  involved  the  adjustment  of  the 
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bolt  locations  to  get  the  best  results.  Clearly,  the  use  of  gap  elements  is  fundamental  in 
obtaining  the  correct  solution.  In  this  respect  die  solution  is  non-linear  since  the  first 
solution  must  define  the  deformation,  and  subsequent  iterations  must  produce  a 
compatible  solution  in  which  the  skin  does  not  'merge'  into  the  spar. 

The  other  area  of  concern  is  at  the  end  of  the  patch,  where  the  patch  picks  up  the  load 
from  the  skin.  Some  bending  and  shear  will  occur  in  this  area. 


5.  Cracked  and  unpatched  model 


5.1  Description 

Firsdy,  the  length  of  the  crack  found  in  aircraft  A8-145  was  48  mm  tip  to  tip  on  the 
inside  of  the  wing,  while  the  length  of  the  crack  on  the  outside  of  the  wing  skin  was 
considerably  shorter  (of  the  order  of  6  mm).  Similar  cracking  has  occurred  to  a  USAF 
F-lllG  in  the  same  area  and  a  fractographic  analysis  has  shown  that  the  crack  had 
started  as  an  elliptical  crack  and  eventually  grew  through  the  outer  surface  of  the  wing 
skin.  In  the  analysis  presented  here  the  crack  will  be  modelled  as  a  through  crack. 
Firstly  this  is  more  convenient,  and  secondly  the  crack  in  the  USAF  F-lllG  was  almost 
through  to  the  outside  surface,  hence  this  assumption  is  only  slightly  conservative. 

The  crack  is  most  conveniently  included  in  the  model  firstly  by  not  applying  restraints 
in  the  x  direction  for  those  nodes  lying  in  the  plane  of  the  crack.  Secondly,  the  elements 
surrounding  the  crack  tip  are  prism  elements  with  the  midpoints  shifted  to  quarter 
points  to  achieve  the  \  /r  stress  singularity. 

5.2  Stress  intensity  factor  results 

The  crack  opening  displacement  method  (COD)  is  used  to  calculate  the  mode  1  stress 
intensity  factor  Ki.  This  is  given  by  ref  (7)  for  plane  strain  as: 

K,  =  (EV/(A. (1  -  v^)))^27i/(La-(L/(2a)))  ...(3) 

and  for  plane  stress  as: 

K,  =  (EV/4.)^27v/(L0.  -  (L /(2a)))  ...(4) 

where  E  =  Young's  modulus 
v  =  Poisson's  ratio 
V  =  displacement 

L  =  length  of  the  element  as  shown  in  Figure  9. 

2a  =  crack  length  (tip  to  tip). 

In  terms  of  fracture  toughness,  the  thickness  of  the  skin  material  corresponds  to  the 
transition  region  of  plane  stress  to  plane  strain.  As  a  result,  the  use  of  the  plane  stress. 
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equation  (3),  is  conservative  because  it  gives  a  higher  value  of  Kj .  Also  for  best  results 
//a  <  1  /20,  however  it  is  known  that  very  small  values  oi  1/ a  may  tmderestimate 
K,. 

The  validation  of  this  formula  has  been  carried  out  using  a  3D  centre  cracked  panel  of 
thickness  of  3.6  mm,  which  is  shown  in  Figure  9.  A  closed  form  solution  is  available, 
ref  (8),  and  for  plane  stress  is  given  by: 

Kj  =  ayfm(l+.l2S(a /b-.2SS(a  /bf  +  \.52?>(a  /bf)  ...(5) 

where  2b  =  the  width  of  the  panel 

2a  =  the  crack  length  tip  to  tip 

This  formulae  is  applicable  in  the  range  0  <a/b  <0.1 
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Figure  9.  Mesh  for  centre  cracked  panel. 

The  F.E.  results  for  a  range  of  mesh  ratio's  denoted  as  'M',  in  Figure  9.  are  given 
Table2.  This  ratio  is  simply  a  means  of  generating  a  gradation  of  elements, 
comparison  of  the  results  has  shown  that  the  error  of  1.4%  has  been  achieved  with  a 
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fine  mesh,  while  a  relatively  coarse  mesh  has  achieved  6%  error.  The  error  is  based  on 
equation  (5)  being  correct.  Also,  the  size  of  the  crack  tip  element  'I'  results  in 
significant  error  for  values  of  l/a<  0.005 .  The  height  of  the  crack  tip  is  defined  as  'h' 
as  shown  in  Figure  9.,  and  the  h/l  ratio  is  also  shown  in  Table  2.  While  meshing 
requirements  for  the  F-111  wing  F.E.  model  have  restricted  the  size  of  the  crack  tip 
element,  the  value  of  /  is  0.3  mm,  and  for  the  48  mm  crack  the  /  /  a  ratio  is  0.0125.  As 
a  result,  an  error  of  less  than  5%  is  expected.  The  dimensions  2a  and  2b  are  500  mm 
and  1000  mm  respectively,  with  a  thickness  of  10  mm. 


Table  2  Validation  of  COD  method  for  determination  of  K, 


Mesh  ratio 

h/l 

1/a 

Kj  (plane  stress) 

(MPa-Jm) 

%  error 

123 

0.5 

0.02 

1.109 

-6.2 

1234 

0.5 

0.02 

1.141 

-3.5 

12345 

0.5 

0.02 

1.186 

-2.2 

123456 

0.5 

0.02 

1.163 

-1.6 

1234567 

0.5 

0.02 

1.165 

-1.4 

12345678 

0.5 

0.02 

1.166 

123456789 

0.5 

0.02 

1.165 

-1.4 

1234 

1.0 

0.01 

1.115 

-3.5 

1234 

1.0 

0.005 

1.052 

-11.1 

The  stress  intensity  formtda  was  therefore  considered  valid,  and  the  crack  elements 
were  introduced  into  the  wing  skin  model. 


Stress  intensities  have  been  evaluated  at  three  points  through  the  wing  skin  thickness. 
These  correspond  to  inner,  middle  and  outside  surface  of  the  skin.  These  values  vary 
slightly,  and  reflect  the  bending  taking  place,  however  only  an  average  value  is 
considered.  Also,  the  result  from  each  crack  tip  varies  sUghtly,  and  in  this  case  it  is 
necessary  to  calculate  the  root  mean  square  stress  intensity  factor.  Figure  24  shows  the 
variation  of  stress  intensity  with  crack  length  for  a  nrunber  of  configurations.  At 
present  only  the  cracked  and  xmpatched  case  (skin  only  and  skin/ spar)  will  be 
considered.  Each  stress  intensity  curve  contains  8  data  points  and  the  cracks  range  in 
length  from  15  -  80  mm,  measured  tip  to  tip. 

Before  the  spar/ skin  model  was  adopted  a  skin  only  model  was  considered.  However 
it  was  clear  that  the  spar,  which  bridges  the  crack,  would  provide  more  restraint  to 
crack  opening  as  the  crack  increased  in  size.  This  effect  is  clearly  seen  in  Fig.  24  in 
which  the  spar  reduces  the  stress  intensity  by  20%  for  large  crack  lengths.  As  expected, 
for  small  cracks,  the  spar  has  less  influence. 

The  value  of  the  stress  intensity  factor  is  dependent  on  the  3D  geometry  and  crack 
length.  In  the  region  of  the  central  stiffener  the  skin/spar  model  gives  values  of 

SOMPayfm  for  crack  lengths  of  15  mm.  As  the  crack  grows  the  stress  intensity 
increases  to  60MPa^/m  at  the  edge  of  the  central  stiffener.  Beyond  this  point  the  stress 
intensity  does  not  change,  presumably  as  a  result  of  the  crack  tip  moving  away  from 
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the  area  of  high  stress.  In  between  tfie  central  stiffener  and  the  auxiliary  stiffener  the 
stress  intensity  rises,  however  it  decreases  to  SQMPayfm  as  it  enters  the  auxiliary 
stiffener.  Once  the  crack  tip  reaches  the  flat  top  surface  of  the  stiffener  the  stress 
inteiisity  increases  again.  Although  cracks  above  80  nun  have  not  been  considered,  it  is 
expected  that  the  stress  intensity  will  rise  rapidly  as  the  crack  grows  outside  the 
auxiliary  stiffener. 

The  value  of  fracture  toughness  given  by  ref  (1)  is  46MPa^fm ,  hence  without  a  repair 
scheme  failure  would  be  predicted  to  occur  at  77%  of  Design  Limit  Load  (DLL). 


6.  Patched  model 


6.1  Description 

The  orientation  of  the  pHes  used  for  the  repair  scheme  is  shown  in  Figure  10.  and  has 
been  obtained  from  ref  (1).  The  patch  consists  of  14  layers  of  Boron  with  a  lay-up  of 
[O2, ±45,03]^ .  The  order  of  the  lay-up  is  given  in  Table  3  where  layer  1  is  the  first  and 

layer  14  is  the  last.  The  patch  has  an  overall  spanwise  dimension  of  468  mm  and  a 
chordwise  dimension  of  324  mm.  Also  the  patch  is  rectangular  but  wifii  cut-off 
comers. 


Figure  10  Lay-up  configuration  for  boron  patch. 
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In  the  F.E.  analysis  the  use  of  symmetry  has  allowed  the  use  of  a  half  model  as  shown 
in  Figure  11.  However,  the  model  still  consists  of  2,000  20  noded  brick  elements  and 
16,000  nodes.  As  a  result  of  computer  size  limitations  and  numerical  accuracy,  there  is 
no  one  model  which  is  used  to  predict  aU  results.  It  has  been  necessary  to  create  a 
number  of  models,  each  with  a  specific  area  of  refinement.  This  has  allowed  a  detailed 
investigation  into  areas  in  which  stress  concentrations  exist. 


Table  3  Ply  lay-up  order 


Layer 

Ply  angle  (deg) 

1 

0 

2 

0 

3 

45 

4 

-45 

5 

0 

6 

0 

7 

0 

8 

0 

9 

0 

10 

0 

11 

-45 

12 

45 

13 

0 

14 

0 

Figure  11.  Spar/skin  F.E.  model  with  outline  of  adhesive  and  boron  patch. 
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6.2  Patch  details 

In  the  F.E.  analysis  the  Boron  patch  is  modelled  as  a  homogenous  orthotropic  material 
in  which  the  material  constants  are  obtained  using  laminate  theory.  The  properties  of 
the  unidirectional  Boron  are  shown  in  Table  4  and  5  while  the  overall  properties  of  the 
laminate  are  shown  in  Table  6,  and  do  not  vary  significantly  over  the  operating 
temperature  of  the  aircraft.  However  compression  strength  properties  do  vary 
sigidficantiy  over  this  temperature  range.  The  material  compliances  for  Boron/  epoxy 
laminate  used  in  the  model  are  shown  in  Table  7.  The  F.E.  outline  of  the  patch  is 
slightly  simplified  as  shown  in  Figure  11.  However  the  F.E.  idealisation  does  include 
the  patch  taper  in  aU  directions.  As  previously  mentioned,  a  'half'  model  is  used  and  is 
shown  in  Figure  11  as  an  exploded  view,  including  the  adhesive. 


Table  4  Mechanical  properties  for  5521  boron/epoxy,  ref  (9) 


Material  property 

Room  temperature 

121®  C 

Tensile  strength  (MPa) 

1520 

1450 

Tensile  modulus  (GPa) 

210 

210 

Compressive  strength  (MPa) 

2930 

1250 

Compressive  modulus  (GPa) 

210 

210 

Flexural  strength  (MPa) 

1790 

1720 

Flexural  modulus  (GPa) 

190 

17 

Interlaminar  shear  (MPa) 

97 

55 

Coeff.  thermal  exp  /®  C 

4.5 

4.5 

Density  g/cm^ 

2.0 

2.0 

Table  5  Properties  for  a  single  layer  of  Boron,  ref  (10) 


=207000.  MPa 

E22  =19000.  MPa 

G,2  =4800.  MPa 

V,2=.21 

a  =  4.1;d0~"  />  C 

Table  6  Overall  properties  of  the  Boron  patch  for  patch  lay-up  o/[02, ±45,03]^ 


=156107. 

EyY=29781. 

=.574884 

=.109675 

=190651. 

Table  7  Compliances  for  Boron/epoxy  laminate  input  to  PAFEC. 


=.640588x10  ^ 

(“  1  /  ^xx) 

=.368264x10  ^ 
i~~^XY  ^ ^Xx) 

5if;^=.524519xl0^(=l/G;^) 

S'yj.  =.335779x10”^ 

(=l/EyY) 

=.705135x10"^ 

^^Yy) 

=.208333x10"^  (=1/G,2) 

S^=.526087xlO~^ 

(=l/^2z) 

=.134524x10”^ 
i~~^xr  ^ ^xx) 

SHyy=.108383x\Q~^ 

(=1/G,2) 
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6.3  Adhesive  details 

The  properties  of  the  FM73  adhesive  are  temperature  dependent  and  are  given  in 
Table  8  for  3  different  temperatures  as  required  for  the  substantiation  of  the  repair.  In 
this  repair  the  adhesive  consisted  of  2  films  with  a  total  thickness  of  0.254  mm  after  the 
repair  application. 


6.4  Modelling  of  Thermal  Effects 

Stresses  and  strains  in  the  various  elements  of  the  repair  are  significantly  influenced  by 
temperature.  This  is  for  two  reasons: 

(a)  The  properties  of  the  adhesive  change  significantly  with  temperature. 

(b)  The  differences  in  thermal  coefficient  of  expansion  for  the  aluminium  compared 
with  the  boron  laminate  results  in  a  thermal  "mis-match"  which  leads  to  residual 
stresses  and  strains. 

The  thermal  mis-match  residual  stresses  can  be  either  beneficial  or  detrimental, 
depending  on  the  location  and  sign  of  the  mis-match.  If  the  coefficient  of  thermal 
expansion  of  the  aluminium  wing  structure  is  greater  than  that  of  the  boron  laminate, 
then  at  ihe  end  of  the  patch  the  thermal  residual  adhesive  shear  stresses  are  opposite  in 
sign  to  ihe  adhesive  shear  stresses  that  arise  as  a  result  of  flight  loads,  and  are 
beneficial.  Over  the  crack  region  however,  the  residual  thermal  adhesive  shear  stresses 
are  the  same  sign  as  the  shear  stresses  due  to  flight  loads,  and  are  detrimental.  If 
however,  the  thermal  coefficient  of  thermal  expansion  of  the  aluminium  wmg 
structure  is  less  than  that  for  the  boron  laminate  then  the  reverse  of  that  just  stated,  is 
true. 

For  the  analysis  presented  in  this  report  the  adhesive  properties  relevant  to  the 
temperature  for  the  case  being  examined  are  always  used.  The  thermal  mis-match  is  a 
significantly  more  complex  analysis  and  it  is  only  included  where  appropriate.  The 
results  detailed  in  section  7  indicate  where  the  thermal  mis-match  effects  have  been 
included  and  why. 
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Shear  stress 
(MPa) 


Shear  strain 


T 


Fig  12a.  Stress  vs  strain  curves  for  FM73,  top;  and  idealised  curve,  bottom. 

The  stress-strain  curves  for  FM73  adhesive  are  shown  in  Figure  12  and  are  taken  from 
the  manufactiures  data  ref  (11).  At  the  lower  temperatures  the  shear  modulus  and 
strength  is  higher  compared  to  that  at  high  temperatures.  However,  at  lower 
temperatures  the  strain  to  failure  is  lower  than  that  at  high  temperatures.  The  failure  of 
an  adhesive  is  governed  by  the  strain  energy  tmder  the  stress  strain  ciurve  taken  to 
failure.  It  is  well  known  that  the  adhesive  transfers  load  primarily  by  shear.  Also,  the 
relatively  low  adhesive  yield  stress  restdts  in  a  considerable  amotmt  of  plasticity. 
However  for  convenience,  the  F.E.  analysis  wiU  be  restricted  to  a  linear  elastic  analysis. 
It  is  possible  to  carry  out  a  strength  analysis  using  a  linear  elastic  analysis  in  which  an 
equivalent  elastic  stress  (or  strain)  is  used  to  predict  failure.  In  this  case  the  stress  (or 
strain)  must  be  based  on  the  strain  energy  under  the  stress  strain  curve  taken  to  fadrure. 
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The  elastic  behaviour  of  an  adhesive  in  a  bonded  joint  usually  involves  a  peak  shear 
stress  at  the  start  of  the  joint,  then  an  exponential  decay.  This  is  the  solution  predicted 
by  the  ID  differential  equation,  ref  (12).  The  F.E.  analysis  however,  shows  a  rapid  rise 
to  a  peak  shear  stress  at  a  distance  approximately  equal  to  the  adhesive  thickness.  In 
the  analysis  presented  here  we  will  examine  the  adhesive  stress  state  at  a  number  of 
critical  points  in  the  patch.  Critical  points  in  the  adhesive  have  been  identified  using 
graphical  means. 


Table  8  Variation  ofFM73  adhesive  properties  with  temperature. 


Temperature 

E(MPa) 

G(MPa) 

V 

-40^  C 

2408. 

892. 

.35 

-32^  C 

2392. 

886. 

.35 

RT 

2273. 

842 

.35 

15°  C 

875. 

324. 

.35 

7.  Results 


7.1  Patched  uncracked 

The  aim  in  this  section  is  to  identify  any  detrimental  effects  that  may  result  from  the 
patching  of  a  structure  which  is  imcracked  or  contains  a  very  small  crack.  Also  it  is 
useful  to  assess  the  ability  of  the  patch  in  being  able  to  stop  the  initiation  of  a  crack. 
The  internal  stresses  and  strains  due  to  the  thermal  mis-match  are  not  evaluated  for 
the  patched  uncracked  model.  The  reason  is  because  the  prime  area  of  concern  is  the 
shear  and  peel  stresses  in  the  adhesive  at  the  end  of  die  patch,  and  Ihe  sign  of  these  is 
usually  opposite  to  the  stresses  which  result  from  the  applied  flight  loads,  i.e.  the 
residual  shear  stress  at  the  end  of  the  patch  is  relieved  by  the  appHed  flight  loads.  In 
terms  of  the  boron  strains,  it  was  useful  to  compare  the  F.E.  strains  with  measured 
strains  from  the  actual  wing  under  load.  The  measured  strains  were  taken  from  a  zero 
set  for  the  cured  patch  in  place,  i.e.  the  residual  strain  was  not  included. 

7.1.1  Boron  strains 

After  the  repair  to  A8-145,  a  strain  survey  was  carried  out  by  the  RAAF,  ref  (13),  in 
which  three  gauges  were  attached  to  the  outer  surface  of  the  boron  patch  to  measure 
spanwise  strains.  The  aircraft  was  loaded  on  the  groimd  by  applying  loads  to  hard 
points  on  the  wing.  One  gauge  was  located  directly  over  the  crack  region  on  the  spar 
centre  line,  while  the  remaining  gauges  were  located  in  the  same  chordwise  plane  liut 
offset  llAmm  either  side  of  the  spar  centre  line.  To  make  valid  comparisons  with  the 
flat  skin  F.E.  model,  the  offset  gauge  results  have  been  averaged.  The  results  have  been 
factored  up  to  the  strain  expected  at  DLL  and  are  shown  on  Figure  13  together  with 
F.E.  predicted  results.  The  centre  gauge  gives  lower  results  than  flie  F.E.  predictions, 
however  the  averaged  results  for  the  offset  gauges,  which  are  in  a  low  spanwise  strain 
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gradient,  agree  to  within  10%.  There  is  however,  some  uncertainty  about  tiie 
comparison  using  these  strain  gauge  results,  since  the  wing  skin  thickness  of  aircraft 
A8-145  is  unknown 

Early  F.E.  models  predicted  strain  levels  of  approximately  7000  microstrain  on  the 
upper  surface  of  the  boron/  epoxy  laminate.  However  the  use  of  two  or  more  elements 
through  the  thickness  of  the  boron  laminate  have  reduced  these  strains  to  levels  of 
2700  microstrain  as  shown  in  Figure  14  and  15.  It  is  however,  not  possible  to  continue 
this  level  of  refinement  to  the  end  of  the  patch.  As  a  result  the  models  with  one 
element  through  the  thickness  of  the  boron  laminate  are  stiU  used  to  predict  boron  and 
adhesives  stresses  where  the  patch  tapers  to  one  layer. 

Another  aspect  of  the  rmcracked  case  is  that  the  boron/ epoxy  patch  picks  up  load  even 
though  a  crack  does  not  exist.  The  strain  level  at  a  distance  60  mm  away  from  the  crack 
is  approximately  2200  microstrain  on  the  lower  surface,  and  it  wiU  be  shown  later  that 
this  does  not  change  even  after  a  crack  is  introduced. 


— □ — Spanwise  strain,  upper 


Spanwise  strain,  lower 


Figure  13.  Chordwise  variation  of  Boron  strain  at  the  crack  location  ,room  temperature  case, 
uncracked  case. 
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— Boron  lower  surfeoe  strain 


.  Boron  strains,  +75'^C,  for  the  uncracked  case. 


—■—Boron  lower  surfeoe  strain  •40°C 
—o— Boron  upper  surface  strain  40°C 


Figure  15.  Spanwise  variation  of  boron  strains,  -4CPCfor  the  uncracked  case. 
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7.1.2  Adhesive  stresses 

7. 1.2. 1  Temperature  of -40  °C 

Adhesive  r^s-  shear  stresses  are  shown  in  Figure  16  and  correspond  to  the  end  region 
of  the  patch.  Figure  16  shows  the  state  of  stress  at  four  different  stations.  The  first 
being  exactly  at  the  end  of  the  patch,  while  the  next  is  0.10  mm  from  the  end  and  so  on. 
The  adhesive  starts  at  a  low  value  and  reaches  a  maximum  value  at  approximately  0.29 
mm  from  the  end  of  the  patch.  Each  curve  has  the  same  characteristic  curve  with  two 
peaks  which  correspond  to  the  chamfered  comers  of  the  patch.  Also  the  stress  value 
corresponding  to  the  location  of  the  forward  auxiliary  spar  (chordwise  co-ordinate  =  0) 
is  considerably  lower  than  the  peak  values.  It  appears  that  the  bending  of  the 
unsupported  plate  away  from  the  spars,  is  also  increasing  the  shear  stress. 

The  Von  Mises  stress  distribution  is  shown  in  Figure  18,  and  has  the  same 
characteristics  as  the  shear  stress  distribution.  The  Von  Mises  stress  distribution  is 
more  representative  of  the  failure  stress  for  the  adhesive  since  it  also  considers  the 
contribution  of  other  stress  components. 

The  peel  stress  distribution  is  shown  in  Figure  17  and  is  a  smoother  distribution  than 
that  given  in  Figure  16  and  18.  However,  the  peel  stress  is  still  higher  where  the  plate 
is  unsupported,  and  clearly  the  bending  of  the  plate  is  responsible  for  this  behaviour. 

Margins  of  safety  wiU  not  be  given  here,  however  in  the  later  sections  the  results 
corresponding  to  cracks,  show  maximum  stress  values  equal  to  those  obtained  in  a 
patched  imcracked  stmcture. 

7. 1.2.2  Temperature  of +75  °C 

The  characteristic  curves  obtained  for  results  of  -40°C  are  also  obtained  at  +75°C  and 
are  shown  in  Figiues  19,  20  and  21.  The  only  difference  being  the  lower  magnitude  of 
the  stresses.  It  will  shown  later,  that  when  a  crack  is  introduced,  the  adhesive  stresses 
at  the  end  of  the  patch  remain  the  same. 
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— A— J\CP  C,  >^^34.95,  uncracked 
— o— -40P  C,  x=234.85,  uncrackad 
— V— Orv  -40^  C,  X=234.66,  uncracked 


Chordwise  coordinate  (mm) 


Figure  16  Shear  stresses  at  end  of  patch  for  -40°  C  case. 


— A— -40^0,  x=234.95,  uncracked 
— ^C,  >^=QM.85,  uncracked 
_^_c^  -40^0,  x:^34.66,  uncracked 


Chordvvise  co-ordinate  (mm) 
Figure  17  Peel  stresses  at  end  of  patch,  ■40°C. 
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V  Mses  (IVRa) 


— A— V  Mses  -40°C,  )f234.95,  uncracked 
— o— V  Mses  -40°C,  x=^.85,  uncracked 
— V— V  Mses  -40°C,  >^^34.66,  uncracked 


-200  --DO  0  IX) 

Chondwise  co-ordinate  (mm) 
Figure  18  Von  Mises  stresses  at  end  of -patch,  -4(PC. 


40 

35 

'sn 


— V— +75^C  )f^34.95 
_a_c^+750C>^^.85 
— V — +7^C  )e234.66 


2  — o— c^+75°C)f^26 


DSTO-TN-0067 


— A— Uncracked,  peel,  +75C)f234.95 
—□—Uncracked,  peel,  +75C)f=234.85 
—■— Uncracked,  peel,  +75Cx=234.66 


Figure  20.  Peel  stresses  at  end  of  patch,  +75'^C . 


— A— Lhcracked  VMses,  +75°C,  >^=234.95 
—o— Uncracked  VMses,  +75°C,  x::234.85 


Figure  21  Von  Mises  stresses  at  end  of  patch,  +75°C. 
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— o — Spanwise  stress  in  Icwe’skin,  patched,  uncracked 
— D— Spanwise  stress  in  upper  skin,  patched,  uncracked 


Figure  22  Spanwise  stresses  in  the  skin  for  the  uncracked  and  patched  case,  DLL,  room 
temperature. 

7.1.3  Stress  distribution  in  the  skin 

The  results  in  Figure  22  show  that  the  stresses  in  the  critical  area  have  been  reduced  to 
approximately  400  MPa,  which  is  the  yield  stress  for  2024-T851.  This  stress  is 
significantly  less  than  the  value  give  in  Figure  8  for  the  unpatched,  uncracked  case, 
and  there  may  be  a  case  for  patching  all  aircraft  in  the  fleet.  Also  in  comparison  to 
Figure  8,  the  backgroxmd  stress  level  (approximately  160  MPa)  under  the  patch  is 
significantly  lower  than  the  applied  stress  of  256  MPa.  Using  graphical  means,  it  has 
been  determined  that  no  significant  detrimental  skin  stresses  have  been  foimd  beyond 
the  end  of  the  patch. 


7.2  Patched  and  cracked 

In  this  section  we  wiQ  consider  the  effectiveness  of  the  patch  and  stress  intensity 
factors  for  a  range  of  crack  lengths.  Also  the  adhesive  and  Boron  stresses  will  be 
considered  for  the  48  mm  crack  length.  The  stress  intensity  factor  work  does  not 
(except  where  indicated)  take  into  accoimt  the  thermal  mis-match  effect,  but  does 
accotmt  for  the  adhesive  properties  at  the  specified  temperature.  A  validation  of  the 
thermal  stress  capability  of  PAFEC  wiU  now  be  carried  out. 
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7.2.1  Validation  of  thermal  stress  analysis 

Experimental  data  has  been  used  to  validate  the  thermal  stress  capability  of  the  PAFEC 
F.E.  program.  In  ref  (14),  an  aluminium  plate  115  mm  x  25  mm  x  1.56  mm  has  been 
reinforced  either  side  with  3  layers  of  a  tmidirectional  Carbon  fibre  and  cured  at  126°C 
using  AS126  adhesive.  Use  of  Carbon  fibre  has  allowed  measurement  of  direct  residual 
stress  using  X-rays.  A  F.E.  analysis  has  been  carried  out  for  this  structure  whose 
material  properties  are  shown  in  Tablel.  The  results  are  shown  in  Figure  23. 
Approximately  6%  error  occturs  between  the  F.E.  and  experimental  values  for  direct 
stress  at  the  midspan  of  the  plate.  The  source  of  the  error  is  not  known.  From  this  data 
the  residual  shear  stresses  can  be  inferred.  Note  the  agreement  of  direct  stresses  at  the 
end  of  the  patch  is  important,  since  this  largely  determines  the  shear  stress 
distribution.  The  experimental  data  points  show  some  scatter  in  the  area  in  which  die 
direct  stress  is  rapidly  changing.  The  plot  shown  in  ref  (14)  does  give  the  95% 
confidence  limit  and  the  F.E.  results  correspond  closely  to  this  plot.  Some  of  the 
experimental  data  on  the  plot  corresponds  to  the  other  end  of  the  structure,  however  it 
is  more  convenient  to  include  this  data  as  shown  in  Figure  23.  The  comparison  with 
the  F.E.  model  and  experimental  data  gives  confidence  for  the  prediction  of  residual 
stress  in  the  lower  skin  repair. 


Fig  23.  Comparison  of  Experimental  residual  data  and  F.E.  results. 
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Table  9  Material  properties  corresponding  to  the  thermal  residual  stress  experiment,  ref  (14). 


Carbon  fibre 

7075-T6 

Adhesive 

Ei=130000.  MPa 

Et  =71016  MPa 

Ga=700.  MPa 

G=7000.  MPa 

Ec  =72397  MPa 

V=0.35 

Ei/E2=15.4 

V=.3 

t=.15 

V=.3 

or,  =0.5xl(H 

a  =23xlCh« 

0. 

7.2.2  Stress  intensity  factors 

Although  tiie  crack  length,  2a,  measured  tip  to  tip,  for  the  repaired  aircraft  was  48  mm, 
a  study  of  crack  lengths  varying  from  15  -  80  nun  was  considered.  At  an  early  stage  in 
the  project  skin  models  were  considered  and  constrained  to  simulate  the  effect  of  the 
spar.  Results  shown  in  Figure  24  show  that  die  skin  models  predicted  high  stress 
intensity  values,  Ki.  The  effect  of  the  spar  is  substantial  in  reducing  the  value  of  Ki, 
particularly  for  large  crack  lengths.  From  Figure  24,  it  is  evident  that  the  RAAF 
designed  patch  is  very  effective  in  reducing  the  stress  intensity  factor  by  almost  50%. 
Using  a  crack  growth  law,  the  reduction  in  crack  growth  is: 

(8a/8N)c/(8a/«N)uc  =  (Kl/Kip)»  ...(1) 

where  for  aluminium  alloys  n  varies  from  3  to  4.  Using  a  value  of  n  =  4  gives  a 
reduction  of  crack  growth  by  a  factor  of  50.  Experimental  results  to  date  have  shown  a 
reduction  in  growlh  by  a  factor  of  100,  however  the  differences  in  the  results  may  be 
due  to  other  experimental  factors.  The  dotted  horizontal  lines  shown  in  Figure  24 
corresponds  to  the  value  of  Ki  predicted  by  C5033,  ref  (1),  for  temperatures  of  -40°C 
and  +75°C,  and  includes  the  contribution  of  thermal  stresses.  Note  the  F.E.  results 
shown  in  Figure  23  do  not  include  thermal  residual  stress. 
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Figure  24  Stress  Intensity  factor  versus  crack  length,  at  DLL. 

The  F.E.  model  was  run  including  the  thermal  effects  for  a  crack  length  of  48  mm  only. 
This  was  done  for  a  range  of  possible  coefficients  of  thermal  expansion  for  the  skin 
material  and  the  results  are  shown  in  Table  10.  This  was  done  at  -40°C  because  this 
condition  gives  die  greatest  temperature  difference  for  the  cure  condition  at  82°C.  At 
an  operating  temperature  of  75°C  the  residual  thermal  stresses  are  neghgible.  From 
Figure  24,  at  the  48  mm  crack  length  and  the  -40°C  condition,  the  F.E.  without  thermal 
effects  gives  a  Ki  of  approximately  13.0MPa4m  .  This  compares  with  a  result  from 
Table  10  with  a  =  15.3jc10”^ /*  C  (as  used  in  C5033)  of  26.7 MPayfin.  The  C5033 
prediction  is  IS.lMPa^fm  .  Including  the  thermal  effect  therefore  increases  the  Ki 
prediction  by  about  16  percent.  The  C5033  prediction  is  considered  to  be  in  good 
agreement  with  the  F.E.  result. 

The  result  in  Table  10  corresponding  to  a  =  23.a:10‘®  /°  C  is  given  and  represents  the 
case  in  which  the  material  value  of  the  coefficient  of  thermal  expansion  of  the  skin  is 
used,  and  over  predicts  the  value  of  Ki.  The  F.E  result  from  Table  10  for 
a  =  6.0;cl0‘'^  /  C  was  included  because  this  is  the  result  which  has  been  derived  from 
strain  gauge  measurements,  ref  (15),  on  a  full  scale  F-111  wing.  The  extreme  values  for 
these  strains  are  shown  in  Table  11.  This  indicates  that  the  value  of  a  assumed  by 
C5033  is  very  conservative,  and  the  residual  thermal  stresses  experienced  in  the  repair 
are  probably  much  lower  than  previously  thought. 
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Table  10.  F.E.  stress  intensity  values  at  -40°C,  corresponding  to  different  coefficients  of 
thermal  expansion. 


a  (/oQ 

23.X106 

15.3x10^(05033  ) 

6.0.xl0^  (Ref  (15)) 

Ki 

MPa  4m 

31.20 

26.68 

22.69 

Table  11  Strain  gauge  readings  corresponding  to  a  simulated  bonding  process. 


Spanwise  strain  near  front  spar 

Spanwise  strain  near  middle  spar 

-100.  |iE 

-370.  |is 

7.2.3  Adhesive  stresses 

The  adhesive  used  is  FM73  which  is  cured  at  a  temperature  of  82°C  and  whose 
properties  are  temperature  dependent.  For  airworthiness  requirements,  a  range  of 
temperatures  must  be  considered  for  the  design  calculations:  -40,  -32  and  75°C.  The 
F.E.  analysis  also  includes  the  contribution  due  to  thermal  stress.  In  the  region  adjacent 
to  the  crack,  the  thermal  residual  shear  stress  in  the  adhesive  is  of  the  same  sign  as  the 
stress  due  to  the  flight  loads  so  the  two  are  added.  In  the  region  at  the  end  of  the  patch, 
the  reverse  is  true  because  the  residual  stress  relieves  the  flight  load.  A  conservative 
approach  is  taken  in  this  report;  the  residual  thermal  effects  are  included  in  the  F.E. 
analysis  in  the  region  of  the  crack,  and  are  ignored  at  the  end  of  the  patch.  The  results 
are  derived  from  different  runs  of  the  same  F.E.  model. 

7. 2. 3. 1  Adhesive  shear  stresses 

Firstly  consider  the  shear  stresses  over  the  crack  region  as  evaluated  by  the  3D  F.E. 
analysis  in  the  mid  plane  of  the  adhesive,  shown  in  Figure  25,  for  the  case  of  -40°C. 
The  boxmdary  condition  at  the  crack  requires  that  X2;x  =  0  at  x  =  0. 


30 


DSTO-TN-0067 


_._F.E  -40PC,  a  f6.24x^0^,  a2=^JCia® 


Figure  25  Adhesive  shear  stress  (z near  crack,  -40°C,  at  DLL. 


Also  it  is  known  tiiat  reaches  a  maximum  at  a  very  small  distance  from  crack.  It 
has  been  necessary  to  have  a  very  fine  mesh,  of  0.05  mm  increments,  to  show  the 
adhesive  shear  starting  at  x  =  0  and  reaching  a  maximum  at  x  =  0.16  mm  of  85  MPa. 
The  procedure  used  in  C5033  is  to  evaluate  the  margins  of  safety  in  the  adhesive,  using 
computations  of  maximrun  strain  at  x  =  0.  For  a  comparison  with  the  F.E.  the  solution 
used  by  C5033,  which  is  based  on  lap  joint  theory  and  assumes  no  bending,  can  be 
formulated  in  terms  of  the  shear  stress  and  is  derived  from  equation  (16)  of  ref  (16): 


r...  =Cl,<T-GJ/(E,IJ  +  a,(GJ/(E,tJ)(a, -a,)hT  ...(6) 

The  first  term  is  due  to  the  applied  stress  and  the  second  term  is  due  to  residual 
thermal  stress. 

Also  the  shear  stress  is  related  to  the  maximum  shear  stress  by: 


^  MAX^ 

These  quantities  are  defined  as: 


...(7) 
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‘=(t. /Gj\E,E,t,t,  /(E,t,  +£.<„>]“ 


...(7) 


where  <t* 


Eq,  tf),  cZq 


g.,k 


AT 

g  jf 


is  the  stress  in  the  absence  of  a  repair 

are  the  boron  modulus,  thickness  and  coefficient  of  thermal 

expansion  respectively 

are  the  skin  modulus,  thickness  and  coefficient  of  thermal 
expansion  respectively 

are  the  shear  modulus  and  thickness  of  the  adhesive 
respectively 

for  a  single  lap  joint  the  thickness  is  divided  by  2. 

is  the  difference  between  the  cure  and  operating  temperatures 

parameters  defined  in  C5033 


The  term  is  a  factor  normally  taken  to  be  1.2  .  This  is  to  cover  the  additional  load 
that  the  repair  attracts  from  the  surrounding  structure  as  a  result  of  the  increased 
stiffness.  T^s  factor  can  also  include  a  bending  correction  term.  However  for 
comparison  purposes  the  lap  joint  results  and  the  F.E.  results  presented  here  are 
tmfactored  and  the  value  of  is  taken  to  be  a  load  attraction  factor  only. 


EquivC5033  75°C.  aj=4.1>dCr®,  02=15.3x10® 


Spanwise  distance  from  cnacK  along  spar  centre  line,  (mm) 


Figure  26  Adhesive  shear  stress  (t^)  near  crack,  +75°C,  at  DLL. 


32 


DSTO-TN-0067 


The  term  Qj.  is  taken  as  a  thermal  inclusion  factor  and  is  taken  to  be  1.5,  see  ref  (1). 
C5033  uses  an  effective  coefficient  of  thermal  expansion  for  the  skin  which  is  a 
measure  of  the  restraint  provided  by  the  surroimding  colder  structure,  and  is  given  by 
=  a(l  +  v)  /I.  The  equivalent  C5033  result  is  plotted  in  Figure  25  and  is  clearly 
conservative  compared  to  the  3D  F.E.  result.  The  +75°C  is  shown  in  Figure  26,  and 
again  the  equivalent  C5033  result  is  conservative.  Also,  in  the  3D  F.E.  result  the 
effective  coefficient  of  thermal  expansion  for  the  skin  material  was  taken  to  be  the 
actual  coefficient  of  thermal  expansion  for  the  material.  This  is  a  conservative 
assumption  in  the  region  of  the  crack  since  the  thermal  contribution  increases  the  total 
shear  stress.  The  maximum  adhesive  stresses  provided  by  the  linear  elastic  analysis 
indicate  that  the  adhesive  is  operating  in  the  plastic  region.  As  previously  mentioned 
in  section  6.3  the  margin  of  safety  may  be  evaluated  using  a  linear  elastic  analysis 
based  on  the  strain  energy  of  the  adhesive  taken  to  failure.  The  equivalent  elastic  shear 
stress  ( Tg )  based  on  diis  strain  energy  is  shown  in  Table  12.  Margins  of  safety  (M.S.) 
have  been  calculated  for  a  range  of  temperatures  and  the  results  are  shown  in  Table  12. 
A  positive  M.S.  exists  throughout  the  temperature  range. 


Tfl&Ze  22  Margins  of  safety  for  adhesive  in  shear  over  crack  for  DUL 


Temp 

G 

Strain  Energy. 

^3Z) 

M.S. 

-40OC 

892 

22.4 

199. 

127.5 

.561 

-320C 

886 

23.6 

205. 

127.5 

.608 

750c 

324 

20.3 

114. 

58.5 

.949 
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— A— ^  -4CPC,  no  residual 

— ^  -4CPC,  x::234.85  no  residual 


Tigure  17.  Variation  of  in  adhesive  at  end  of  patch,  -4(PC. 


Now  consider  the  adhesive  stress  state  at  the  end  of  the  patch,  conservatively  ignoring 
assumption  made  is  that  the  contribution  due  to  thermal  residual  stresses  since  they 
relieve  the  stresses  due  to  flight  loads  at  the  end  of  the  patch.  Figure  27  shows  the 
chordwise  variation  of  the  midplane  stresses,  and  corresponds  to  the  adhesive 
material  property  state  at  -40°C.  Figure  27  is  comprised  of  four  curves  each  at  an 
increasing  distance  from  the  end  of  the  patch.  The  curve  having  the  highest  value  is 
located  approximately  0.30  mm  from  the  end  of  the  patch.  The  peak  value  on  each 
curve  corresponds  to  the  chamfered  comer  of  the  patch,  while  in  comparison,  the 
values  corresponding  to  the  forward  auxiliary  spar  location  (y  =  0)  are  significantly 
lower.  Furthermore,  the  peak  values  are  superimposed  on  a  curve  which  is  itself 
higher  than  those  values  at  the  spar  location.  This  suggests  that  bending  of  the 
unsupported  skin,  midway  between  the  forward  auxiliary  spar  and  centre  or  front 
spar,  is  increasing  the  adhesive  shear  stress  significantly. 

In  the  case  of  a  3D  adhesive  there  are  components  other  than  t^x  that  must  be 
considered.  The  Von  Mises  stress  is  a  better  representation  of  the  failure  stress  for  the 
adhesive,  and  this  distribution  is  shown  in  Figure  28.  Again  the  peak  values  occur  at 
the  chamfered  comers  while  lower  values  exist  in  the  region  of  the  forward  auxiliary 
spar.  Also,  the  maximum  values  occur  approximately  0.3  mm  from  the  end  of  the 
patch. 
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The  corresponding  curves  for  the  adhesive  at  +75°C  are  shown  in  Figures  29  and  30, 
and  are  similar  to  the  stress  plots  corresponding  to  -40°C.  However  for  the  +75°C  case 
the  stress  levels  are  considerably  reduced  as  a  result  of  the  lower  shear  modulus,  and 
again  tire  maximum  values  occur  0.3  mm  from  the  end  of  the  patch.  Results  for  the  - 
32°C  case  are  not  shown,  however  the  shear  modulus  for  the  adhesive  at  this 
temperature  is  almost  identical  to  the  -40°C  case. 


From  Figure  12a  it  is  evident  that  over  the  complete  temperature  range  the  adhesive 
will  not  be  operating  in  the  elastic  regime,  particularly  in  the  region  of  the  chamfered 
comers.  However,  the  margins  of  safety  have  been  computed  for  operation  at  these 
temperatures  and  are  shown  in  Table  13.  Positive  margins  of  safety  exist  throughout 
the  temperature  range.  The  shear  stress  predicted  by  the  ID  lap-joint  equation  are 
shown  in  Table  14.  These  values  are  considerably  lower  than  those  predicted  by  the  3D 
F.E.  analysis  shown  in  Figures  27  &  29. 


Table  13  Margins  of  safety  for  adhesive  at  end  of -patch  for  DUL 


Temp 

G 

S.S. 

T  =  o-y^f  / 1.732 

M.S. 

-40OC 

892 

22.4 

199. 

88.3 

1.25 

-320C 

886 

23.6 

205. 

88.3 

1.32 

750c 

324 

20.3 

114. 

59.5 

0.92 

Table  14  Shear  stress  predicted  by  lap-joint  equation  (MPa) 


-400C 

+750C 

27.48 

16.56 
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— A— VI\/1-40^C,  no  residual,  x=234.95 
—a—\JM-4CPC,  no  residual,  >^234.85 
— V— VM-40^C,  no  residual,  x=234.66 


Chordwise  oo-ordinaie  (mm) 
Variation  of  <jy^  in  adhesive  at  end  of  patch,  -4(PC. 


— A— SZS  +75PC,  x=234.95  resid 
— SZS  +75°C,  >?234.85  resid 
— V— SZE  +75PC,  )?234.66  resid 


Chcxdwise  cxKDitlinate  (mm) 

Figure  29.  Variation  of  r^x  in  adhesive  at  end  of  patch,  for  75°C. 
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_^_Wl+75°Cpf^.95 
— \4VI+75°Cpf234.85 


—A— Peel  -40^0,  >^=234.95  no  residual 
— D— Reel  -4CPC,  >e^!34.85  no  residual 
—'7— Pee!  -40^C,  >f=234.66  no  residual 


Figure  31  Variation  of  a  in  adhesive  at  end  of  patch,  for  -4(FC. 
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— A— Pgel  stress  +75PC,  x=234.95 
—□—Peel  stress  +75°C,  >^^34.85 


7.23.2  Adhesive  peel  stresses 

The  acJhesive  peel  distribution  around  the  end  of  the  patch  is  shown  in  Figure  31  for 
tiie  case  of  ■40°C.  Peak  values  occur  in  the  region  of  the  chamfered  comer  while  lower 
values  occur  in  the  region  of  the  spar  (y  =  0).  The  maximum  peel  stress  occurs  exactly 
at  the  end  of  the  patch  and  decays  rapidly  in  a  distance  of  0.7  mm.  Also  the  slightly 
higher  values  away  from  the  spar  location  indicate  that  bending  of  the  unsupported 
plate  has  a  slight  influence.  The  peel  stress  distribution  corresponding  to  adhesive 
material  properties  at  +75°C  are  shown  in  Figure  32.  In  comparison  to  the  -40°C  case 
the  stresses  are  lower  simply  due  to  the  lower  adhesive  shear  modulus  at  the  higher 
temperature. 

Comparisons  of  the  F.E.  are  made  with  the  d.e.  expression,  given  by  Hart-Smith,  ref 
(12).  This  expression  is  used  by  C5033  and  can  be  used  for  a  single  lap  joint  but  with  no 
bending: 


Op  =  -Wo  Ec/ta  ...(9) 

where 

Wo  =  (A/ 2)  (e  X(x-l/2)  +  e  -X(x  +  1/2))  sin  Xx 

+(B/2)  (e  X(x-l/2) .  e  -X(x+  1/2))  Cos  lx  ...(10) 

where 
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X  =  (Ec/(4Dta))0-25 
A  =  (x(ta/2D)SinX//2)/x3 
B  =  -(x  (ta/2D)  Cos  X  //2)/  X^ 

=  effective  tensile  modulus  of  the  adhesive 
X  =  shear  stress  (normally  taken  to  be  the  plastic  value  tp) 

D  =  Eq  to^/12(l-v2)  and  is  the  bending  stiffness  of  boron 

In  this  case  the  location  of  the  coordinate  system  is  at  the  midpoint  of  the  joint,  and  x  = 
±1/2  corresponds  to  the  ends  of  the  joint. 

A  comparison  of  the  d.e.  results  and  F.E.  work  is  shown  in  Figure  33  and  34.  Again  the 
difference  between  the  results  is  substantial.  Secondary  bending  results  in  very  large 
peel  stresses  which  exceed  the  allowable  strength  of  70mpa,  ref(l),  over  a  very  small 
region  of  less  than  0.5  mm  diameter.  While  local  failure  of  the  adhesive  may  occur  this 
does  not  mean  that  in  this  particular  case  that  the  patch  itself  wiU  fail.  A  disbond  of  6.0 
mm  was  considered,  and  as  a  resxilt  the  peel  stress  had  reduced  to  that  predicted  by 
the  d.e. 

In  deriving  a  solution  for  the  d.e.  for  the  peel  it  is  normally  assumed  that  the  shear 
stress  is  constant.  Generally  this  is  the  case  when  looking  at  the  maximum  strength  of  a 
bonded  joint,  where  the  adhesive  is  in  the  plastic  region. 


Equiv.  C5033  peel  stress  crack  region 
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—■—RE  pee!  stress  75°C,  crack  region 


Figure  34  Spanwise  variation  of  peel  stress  away  from  crack. 


— A— VMses  square  234.95  -40°C 
— □— VMses  square  234.85 
—V— VMses  square  234.66  40°C 


Figure  35  Von  Mises  stresses  at  the  end  of  the  patch  for  square  comer  configuration, -4(PC. 
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—A— Peel  square  cnr  234.95 
— D— Pee!  square  cnr 234.85 
— V— Peel  square  cnr 234.66  -40°C 


Figure  36  Variation  of  at  the  end  of  a  patch  with  a  square  comer  configuration,  for 

temperature  of-4CPC. 


7. 2.3.3  Comparison  of  chamfered  and  square  comer  configuration  patch 

A  comparison  has  been  carried  out  to  determine  the  advantage  of  using  a  square 
patch.  Figure  35  shows  the  Von  Mises  stress  distribution  in  the  adhesive 
corresponding  to  a  square  comer  patch  at  -40°C,  and  is  directly  comparable  to  Figure 
28  for  the  chamfered  comer  configuration.  Given  that  for  a  3D  stress  state  the  Von 
Mises  stress  is  the  best  representation  of  failme  for  an  adhesive,  this  indicates  that  the 
square  comer  patches  may  have  advantages.  Ftirihermore,  it  has  been  foimd  in  terms 
of  stress  concentration  in  the  skin,  the  square  comer  configuration  resxilts  in 
significantly  lower  stress  values.  In  terms  of  peel  stress,  a  comparison  of  Figures  36 
and  Figure  31  shows  no  discernible  difference. 

7.2.4  Boron  strains 

As  a  result  of  the  low  effective  coefficient  of  thermal  expansion  in  the  wing  skin, 
thermal  residual  stresses  in  the  boron  have  not  been  included.  Figure  37  shows  the 
chordwise  variation  of  the  spanwise  boron  strains  either  side  of  the  spar  centre  line  for 
the  lower  surfaces.  The  dip  in  the  strain  distribution  at  the  spar  centre  line  indicates 
that  secondary  bending  is  present.  These  curves  are  similar  to  those  shown  for  the 
uncracked  case  shown  in  Figure  13.  The  behaviour  of  the  patch  directly  over  the  crack 
area  is  clearly  sensitive  to  temperature.  The  minimum  strains,  directly  over  the  spar  for 
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tihe  ■40°C  and  room  temperature  case  are  much  lower  than  the  75°C  case.  An 
explanation  is  that  at  75°C  the  lower  stiffness  of  the  adhesive  is  unable  to  constrain  the 
boron  and  less  bending  occurs  at  this  point  The  measured  boron  strains,  ref  (13),  have 
been  scaled  up  to  DLL  and  are  shown  in  Figure  37.  The  F.E.  predicted  strains  away 
from  the  spar  centre  line  are  10%  higher  tiian  the  RAAF  measrires  strains,  however  the 
actual  skin  thickness  of  the  aircraft  involved  is  not  known.  To  evaluate  the  wing  skin 
thickness  effect  on  these  strains  a  wing  skin  thickness  of  4.0  mm  was  considered  as 
shown  in  Figure  38.  The  agreement  widi  measured  strains,  located  away  from  the 
centre  line,  has  improved  significantly.  The  centre  gauge  (gauge  length  of  12.5  mm)  is 
in  a  region  of  a  very  high  strain  and  when  the  F.E.  results  are  averaged  over  the  length 
of  the  strain  gauge  some  error  exists.  Ultrasonic  measurements  of  the  thickness 
distribution  in  the  rmcracked,  tmpatched  wing  indicate  that  significantly  more  skin 
material  has  been  removed  in  the  spanwise  direction  than  is  allowed  for  in  the 
drawings.  As  a  result  the  spanwise  strain  gradient  would  be  reduced.  If  this  is  also  the 
case  with  the  repaired  aircraft  then  this  would  explain  the  error  obtained  at  the  centre 
gauge. 

The  results  shown  in  Figure  39  correspond  to  spanwise  strains  in  the  boron  along  the 
spar  centre  line  for  the  lower  surfaces  and  extend  to  a  distance  of  85  mm  from  the 
crack.  These  strains  are  approximately  3000  microstrain  and  reach  a  peak  at  10  mm 
outboard  of  the  crack  along  the  spar  centre  line.  Results  shown  in  Figure  40  correspond 
to  the  upper  surface  strains  whose  peak  values  are  approximately  11,000  microstrain 
directly  over  the  crack,  at  DLL.  In  this  region  two  elements  have  been  used  through 
the  thickness  of  the  boron,  and  as  a  result  the  predicted  strains  are  considered  to  be 
acctirate.  These  strains  are  considerably  higher  than  those  obtained  in  the  uncracked 
case  shown  in  Figure  15.  It  is  likely  that  this  increase  is  due  to  the  change  of  the  load 
path  in  the  region  of  the  crack.  The  load  is  now  carried  by  the  boron  and  the  neutral 
axis  offset  has  been  increased,  resulting  in  increased  bending  of  the  boron. 

It  has  been  shown  previously  that  the  peel  stresses  in  the  adhesive  in  the  region  of  the 
crack  will  result  in  local  failure.  A  representative  test  article,  used  in  the  structural 
testing  program,  has  shown  disbonds  of  approximately  5  mm  either  side  of  the  crack 
vmder  representative  loads.  A  6  mm  disbond  has  been  included  in  a  F.E.  model,  and 
the  corresponding  boron  strains  are  shown  in  Figure  41.  The  highest  strains  occur  5 
mm  away  from  the  crack  and  are  approximately  4500  microstrain.  It  is  likely  that  the 
plasticity  of  the  adhesive  at  the  disbond  will  reduce  these  strains  and  the  maximum 
values  considered  are  those  occurring  over  the  crack.  The  disbond  has  also  reduced  the 
magnitude  of  the  lower  surface  strains  for  both  -40°C  and  75°C  temperatures. 

Maximum  strains  in  the  boron  for  the  required  temperature  range  and  margins  of 
safety  are  shown  in  Table  13. 
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— A— ^Danwise  boron  strain  +75°C 
—+—^Danwise  boron  strain  RT 


Figure  37  Spanwise  lower  surface  boron  strain  in  plane  of  crack  for  ■40°C,  room  temperature 
and  75°C  cases. 


Boron  spanwise  strain,  lov\er  t=4.0mm 
■  FW^nneasuiad  strains  e>4iap.  DLL 


Figure  38  Spanwise  lower  surface  strain  in  plane  of  crack  for  -4(PC  corresponding  to  a  skin 
thickness  of  4.0  mm. 
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—■—Boron  spanwise  strain  lov\er  +75°C 

0.0035 

0.0030 

0.0025 

0.0020 

0.005 

0.005 


Spanwise  distance  trom  crack  (mm) 
Spanwise  boron  lower  surface  strain  for  48mm  crack. 


Spanwise  distance  from  crack  (mm) 
Figure  40  Spanwise  boron  upper  surface  strain  for  48  mm  crack. 
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—■—disbond,  boron  spanvuise  strain,  upper  75°C 
—o— disbond,  boron  spanwise  strain,  upper  •40°C 
—•—disbond,  boron  spanwise  strain,  lower  75°C 


Figure  41  Disbond,  variation  of  spanwise  upper  surface  boron  strain  with  distance  from  crack. 


Table  15.  Margins  of  safety  in  boron  at  DUL. 


Temperature  ®C 

Maximum  allowable 
strain 

Ultimate  strain 

m.s. 

-40 

6000. 

6150. 

-.024 

-32 

6000. 

6150. 

-.024 

+75 

6000. 

6750. 

-.111 

7.2.5  Boron  interlaminar  shear 

Graphically,  the  maximum  interlaminar  shear  stresses  have  been  fotmd  to  occur  over 
the  crack.  The  allowables  are  given  in  Tablel6,  together  with  the  maximum  values  in 
the  repair  at  DUL.  The  margins  of  safety  are  positive  throughout  the  temperature 
range. 


Table  16  Interlaminar  shear  in  region  of  crack 


Temperature  “C 

Allowable  Int  shear  (MPa) 

Interlaminar  shear  (MPa) 

m.s. 

-40 

97 

41.16 

1.36 

-32 

97 

41.16 

1.36 

+75 

55 

22.2 

1.48 
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8.  Comparison  of  C5033  and  3D  F.E. 


8.1  Source  of  equations  used  in  RAAF  design 


The  following  formulae  used  for  patch  design,  have  been  compared  with  the  results  of 
the  3D  F.E.  .  While  the  derivation  of  these  equations  will  not  be  carried  out,  where 
possible  the  source  of  these  equations  will  be  identified. 

8.1.1  Stress  in  the  skin  under  the  repair 

Following  is  the  equation  for  the  direct  stress  in  the  skin  tmder  the  repair: 

o-„  =€l,^(7/(\  +  s)  +  Q.TE^s/(\.  +  s))(a^  -a2)\Sl\  -(11) 

where  s  =  / {E^t^)  —(12) 

=  1.44 

=  1.5 

This  equation  is  derived  from  equilibrium  considerations  and  is  made  up  of  2  terms. 
The  first  term  identifies  the  stress  level  in  the  skin  resulting  from  load  transfer  into  the 
region,  while  the  second  term  is  a  result  of  the  thermal  stress.  Normally  the  load 
attraction  factor  =  1.2 .  In  this  case  however,  two  correction  factors  have  been 
applied  to  .  One  is  due  to  the  increased  load  the  patched  structure  attracts,  due  to 
increased  stiffness.  The  other  is  caused  by  a  stress  concentration  due  to  bending 
effects,  the  derivation  of  which  is  included  in  ref  (1).  Consider  now  the  second  term  of 
equation  (11).  This  is  the  direct  stress  in  the  skin  as  a  result  of  thermal  stress/residual 
stress,  and  is  derived  from  equation  (6.1a)  of  ref  (17).  Also,  this  term  contains  a 
correction  factor  of  Qj.  =  1.5 , .  Thermal  stresses  may  be  considered  to  be  bi-directional 
in  the  plane  of  the  adhesive,  and  hence  some  stress  concentration  factor  is  appropriate 
for  bending  effects. 

8.1.2  Maximum  adhesive  shear  strain 

The  expression  for  maximum  adhesive  shear  strain  is  based  on  simple  lap  joint  theory 
which  has  been  extended  to  include  plastic  deformation.  The  first  step  used  in  ref  (1) 
was  to  compute  the  shear  strain  which  defines  the  elastic  limit.  This  equation  is 
derived  from  equation  (5.26)  of  ref  (16),  and  is  given  by: 


Ymax=^M/Ga 


...(13) 


If  this  exceeds  the  elastic  limit  for  the  adhesive  shear  strain,  then  the  plastic  limit  must 
be  computed  and  is  derived  from  equation  (5.27)  of  ref  (16)  and  is  given  as: 
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r.„=(f,/(2G.)|l  +  (<T.»,A/r,f]  -(14) 

where  the  parameter  A ,  is  given  by  equations  (5.9)  of  ref  (16)  as; 

1 = [(G.  /«4i  /(£,o + 1  /(Vo)]p  -(15) 

where  the  elastic  transfer  length  is  given  by  1  /  A . 

Essentially,  the  margin  of  safety  is  based  on  the  computation  of  available  energy 
calculated  as  the  area  xmder  the  adhesive  stress  strain  cxirve  taken  to  failure. 

8.1.3  Stress  intensity  in  the  repaired  skin 

The  stress  intensity  in  the  repaired  structure  is  a  function  of  the  stress  in  the  skin  tmder 
the  patch,  and  other  patch  parameters.  The  concept  of  an  upper  boxmd  for  ,  which 
is  independent  of  crack  length,  was  derived  from  equations  (5.22)  and  (5.28)  developed 


in  ref  (16)  and  is  given  by: 

...(16) 

where  =1  +  2(a Qt^X !  t p) 

...(17) 

and  4=1  +  (o-^iiA  / 

...(18) 

The  results  shown  in  Figure  13  which  correspond  to  the  skin/  spar  patched  results, 
indicates  that  the  value  of  is  constant  and  independent  of  crack  length. 

8.1.4  Adhesive  peel  stresses 

The  derivation  of  the  maximtun  peel  stress  is  given  by  equation  (83)  of  ref  (12)  as  a 
fimction  of  the  shear  stress: 

=  r,[3E‘^t,a-y")/(EA>r’  -(W) 

where  E'^  is  the  effective  transverse  tensile  modulus  of  die  adhesive  and  is  given  by 
equation  (68)  of  ref  (12)  as: 

4  =  (1  /Ec  +2/E,+4 /E,y  ...(20) 
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The  value  of  Xp  is  taken  to  be  tiie  ultimate  failure  strength  of  the  adhesive.  The 
calculation  of  the  peel  stresses  have  been  evaluated  in  the  region  of  the  crack  and  at 
the  end  of  the  patch. 

8.1.5  Maximum  stress  in  patch 

The  maximum  direct  stress  is  derived  from  simple  equilibrium  considerations  in 
which  no  bending  occurs  and  is  given  by: 

<7^=atJ  -(21) 

This  formulae  is  only  applicable  to  the  region  directly  over  the  crack.  Computations  of 
interlaminar  shear  for  Boron/ epoxy  are  not  considered. 

8.1.6  Strength  properties 

A  summary  of  ultimate  strength  properties  obtained  from  ref  (1)  is  shown  in  Table  8. 
The  M.S.  derived  from  F.E.  results  are  based  on  these  strengths. 


Tahle  17  Strength  properties  of  skin,  boron/epoxy  and  adhesive. 


Temperature 

Temperature 

Temperature 

Parameter 

-400C 

-32°C 

750c 

Skin  ultimate  strength  in  the 
rolling  direction  2024r-T851  (MPa) 

448 

448 

448 

Skin  fracture  toughness 

Kjc  (MPayfm ) 

46.2 

46.2 

46.2 

Adhesive  shear  strain  Ye 

.0903 

.0866 

.067 

Adhesive  strain  shear  y  p 

.4377 

.486 

1.493 

Adhesive  peel  stress  (MPa) 

70.0 

70.0 

70.0 

Boron  laminate  ultimate  strain  / 

.006 

.006 

.006 

The  values  shown  in  Table  17  are  derived  from  the  RAAF  design  package,  ref(l),  and 
rely  on: 

1.  2024-T851  Properties  from  ref  (5) 

2.  Fracture  toughness  tests  at  AMRL 

3.  FM73  Manufacturer's  data,  ref  (11) 

4.  Boron  Manufactures  data,  ref  (9) 

The  peel  strength  allowable  of  70  MPa,  ref(l),  used  throughout  the  temperature  range, 
is  open  to  question.  Data  for  FM73  at  room  temperature  (RT),  from  ref(19),  gives  an 
ultimate  tensile  strength  of  44.2  MPa  and  is  taken  to  be  the  ultimate  peel  strength. 
Unfortimately  data  does  exist  for  other  temperatures.  Simple  ratios  have  been  fotmd  to 
exist  between  shear  yield  and  tensile  yield  for  epoxy  adhesives  in  refs  (20)  and  (21). 
Extending  this  to  ultimate  strengths  and  using  data  from  Figure  12  and  42,  some 
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approximate  peel  strengths  have  been  computed  as  shown  in  Table  18.  (The  shear 
strengths  for  -40‘’C  and  -32°C  are  assumed  to  be  approximately  the  same  as  the  -550C 
case).  These  peel  strengths  are  based  on  a  ratio  of  peel/ shear  ratio  of  44.2/40  and  are 
derived  from  the  RT  data.  While  the  RT  ultimate  strength  is  44.2  MPa,  the  equivalent 
elastic  failure  strength  based  on  the  same  strain  energy  to  failuxe  is  55.4  MPa,  as  shown 
in  figure  42.  A  further  factor  of  (55.5/44.2)  is  applied  to  the  ultimate  peel  strengths  to 
obtain  the  equivalent  elastic  allowable,  and  is  shown  in  Table  18.  There  is  clearly,  a 
need  for  temperature  dependent  peel  stress  data  for  FM73. 


Table  18.  Estimated  peel  strength  variation  with  temperature  (also  includes  RT  data) 


Temperature  °C 

-40 

RT 

+75 

Estimated  peel  strength  (MPa) 

63.5 

44.2 

37.2 

Eqmvalent  Elastic  allowable  (MPa) 

81.6 

55.5 

46.8 

- R\/73  strain  rate  10^/sec 

—+— Equivalent  linear  elastic 


Figure  42  Tensile  stress  strain  curve  for  FM73  at  room  temperature  and  at  a  strain  rate  of 
10^/sec,  from  ref (19). 


8.2  Summary  of  3D  F.E.  stress/strain  results 


A  summary  of  the  3D  F.E.  results  are  shown  in  Table  18. 

8.2.1  Skin 

The  F.E.  analysis  confirmed  that  the  repaired  structure  containing  a  48  mm  crack  has  a 
positive  M.S.  at  DUL,  based  on  a  fracture  toughness  value  (Kj^.)  of  46  MPaVm. 
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Furthermore  the  associated  testing  program  has  confirmed  the  value  of  Kj(^  .  The 
margin  of  safety  at  +75°C  is  slightly  higher  than  at  the  lower  temperatures,  mainly  as  a 
result  of  the  lower  thermal  residual  stress.  The  analysis  has  also  confirmed  that  a 
positive  M.S.  exists  at  DUL  for  the  stress  in  the  skin  throughout  die  operating 
temperature  range. 

8.2.2  Adhesive 

8.2.2. 1  Adhesive  shear 

The  maximum  shear  stresses  for  the  -40°C  case  have  been  foimd  to  be  near  the  crack 
region,  while  for  the  +75°C  case  die  stresses  at  the  chamfered  comer  are  equal  to  those 
near  the  crack  region. 

8. 2.2.2  Adhesive  peel 

The  maximum  peel  stresses  occur  at  the  crack,  and  are  caused  by  large  out  of  plane 
displacements.  While  local  failure  will  occur  near  the  crack  over  a  very  small  region, 
this  failxue  will  not  spread.  The  important  peel  stresses  are  those  at  the  end  of  the 
patch  since  failtue  in  this  region  wiQ  result  in  complete  failure  of  the  patch.  M.S.  are 
based  on  approximate  data  and  failure  at  the  lower  temperatures  is  close  to  the  DUL. 

8.2.3  Boron 

8. 2. 3.1  Maximum  strains 

The  maximmn  strains  occur  over  the  crack  region  of  the  upper  sruface  of  the  boron, 
and  is  equal  to  12000  microstrain  at  DLL.  However  a  disbond  is  expected  to  form  in 
this  area  before  DLL  is  reached.  The  maximum  expected  strain  in  the  boron  at  DUL, 
after  a  disbond,  is  6750  microstrain.  The  maximum  lower  surface  strains  occur  a 
distance  of  10  mm  away  from  the  crack  is  equal  to  approximately  3800  microstrain.  A 
small  negative  M.S.  occurs  throughout  the  temperature  range. 

8.2.3. 2  Interlaminar  shear 

The  maximum  interlaminar  shear  stress  occurs  over  the  crack  region  and  a  positive 
M.S.  exists  throughout  the  operating  temperature  range. 
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Table  18.  A  summary  of  the  3D  F.E.  results  at  DUL 


Parameter 

-40OC 

-320C 

750c 

Skin:  Stress  intensity 
factor  (patched) 

(  MPa-Jtn ) 

26.68x1.5=40.0 

26.68x1.5=40.0 

25.10x1.5=37.7 

Adhesive:  Shear  over 
crack.(MPa)(Table  12) 

1.5x118.  =  153. 

1.5x118.  =  153. 

1.5x67.=  103. 

Adhesive:  Peel  over 
crack(MPa) 

230x1.5  =  345. 

230x1.5  =  345. 

87x1.5  =  131. 

Adhesive:  Von  Mises  at 
end  of  patch. 

(MPa)(Tablel3) 

1.5x1.732x88.3  =  229.4 

1.5x1.732x88.3=  229.4 

1.5x1.732x59.5=154.6 

Adhesive:  Peel  at  end 
of  patch  (Fig.  30,31) 

59.3x1.5=88.95 

59.3x1.5=88.95 

29.2x1.5=43.8 

Boron:  Direct  strain 
over  crack  (Table  15) 

1.5x4100  =  6150.U 

1.5x4100  =  6150u 

1.5x.4500=6750  .u 

Boron:  direct  strain  at 
end  of  patch 

1.5x1730  =  2595u 

1.5x1730=  2595u 

1.5x2930=  4395u 

Boron:  Inter-laminar 

shear,  stress  adjacent  to 
crack  (Table  16) 

approx  1.5x27.44  = 

41.16  MPa 

41.16  MPa 

approx  1.5x14.8  =  22.20 
MPa 

Boron:  Inter-lam.  shear 
stress  at  end  of  patch 

approx  1.5x12.  =  18.0 
MPa 

approx  1.5x12.  =  18.0 
MPa 

approx  1.5x8.89  =  13.33 
MPa 

8.3  Comparisons  of  margins  of  safety 

The  overall  results  are  shown  in  Table  18.  A  comparison  of  the  results  from  the  F.E. 
analysis  and  results  predicted  by  C5033,  are  shown  in  Table  19.  Note  that  the  M.S. 
originally  calculated  by  the  RAAF  in  ref  (1)  correspond  to  DLL.  The  following  values 
have  been  re-calctdated  and  now  correspond  to  DUL. 

8.3.1  Stress  intensity  factor 

A  small  difference  is  noted  between  the  M.S.  for  stress  intensity  factor  predicted  by 
F.E.  and  C5033.  The  F.E.  results  show  a  lower  M.S.  for  Kj  at  the  lower  temperatures 
and  is  due  to  lower  residual  tiiermal  stress  experienced.  However,  at  a  temperature  of 
+75°C  die  M.S.  is  higher. 

8.3.2  Adhesive  shear 
8.3.2. 1  Over  crack 

The  M.S.  for  adhesive  shear  over  the  crack  predicted  by  the  F.E.  method  is  higher,  and 
is  largely  due  to  the  conservative  value  used  for  Qj  in  ref  (1). 


8. 3. 2. 2  End  of  patch 

The  M.S.  for  adhesive  shear  at  the  end  of  the  patch  is  not  computed  in  ref  (1),  since  it  is 
assumed  that  the  critical  shear  values  occur  in  the  region  of  the  crack. 
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8.3.3  Adhesive  peel 

The  M.S.  for  adhesive  peel  predicted  by  the  F.E.  method  at  the  end  of  the  patch  is 
significantly  different,  since  peel  stress  failure  data  used  is  dependent  on  temperature. 

8.3.4  Boron  strain 

The  M.S.  predicted  by  the  F.E.  for  the  boron  strain  is  substantially  lower  and  is  due  to 
bending  in  the  region  of  the  crack,  as  a  result  of  the  neutral  axis  offset.  The  F.E.  is 
conservative  because  it  did  not  take  into  accotmt  the  thermal  residual  compression 
strain  which  is  beneficial. 

8.3.5  Overall  M.S. 

While  some  disagreement  occurs  over  the  values  of  M.S.  it  is  clear  that  the  structural 
integrity  of  the  RAAF  designed  repair  is  confirmed. 


Table  19.  Comparison  ofm.s.  predicted  by  C5033  and  3D  F.E.  for  48  mm  crack  at  DUE. 


C5033 

C5033 

C5033 

3D  F.E. 

3DF.E. 

3D  F.E. 

Parameter 

-400  C 

-320C 

750  c 

-400  C 

-320C 

750  c 

Skin:  Stress 

intensity  factor 
(patched) 

0.227 

0.087 

0.019 

.155 

.155 

.225 

Adhesive:  Shear 
over  crack. 

Table  10 

0.448 

0.399 

1.581 

0.561 

0.608 

0.949 

Adhesive:  Peel 
over  crack 

- 

- 

- 

* 

* 

* 

Adhesive:  shear 
at  end  of  patch.  . 
(Von  Mises) 

Table  11 

1.25 

1.25 

0.92 

Adhesive:  peel 
at  end  of  patch 

0.349 

0.390 

5.271 

-.083 

-.083 

.068 

Boron:  direct 
strain  over 

crack 

0.610 

0.359 

0.359 

-0.024 

-0.024 

-0.111 

Boron:  direct 

strain  at  end  of 
patch 

■ 

0.092 

0.093 

0.128 

Boron:  Inter- 
lam.  shear 

stress  adjacent 
to  crack 

1.36 

1.36 

1.48 

Boron:  Inter- 
lam.  Shear 

stress  at  end  of 
patch 

4.39 

4.39 

3.13 

*  failure  limited  to  a  small  region 
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9.  Conclusions 


A  3D  F.E.  model  of  the  skin/ spar  structure  has  been  developed  and  validated  with 
strain  gauge  data  from  a  fuU  size  tmcracked  wing.  The  3D  F.E.  model  has  been 
extended  to  model  the  fatigue  crack,  the  RAAF  designed  boron  fibre  patch  and  refined 
to  the  extent  necessary  to  evaluate  adhesive  parameters.  Furthermore  strain  gauge 
readings  made  by  the  RAAF  on  the  repaired  wing  have  also  validated  the  F.E.  model 
with  the  composite  repair. 

Stress  intensity  factors  have  been  evaluated  for  a  range  of  crack  lengths  for  the 
imrepaired  and  repaired  structure,  and  the  crack  growdi  rate  is  expected  to  be  reduced 
by  a  factor  of  100,  and  is  in  agreement  with  test  results. 

The  reduction  of  the  maximum  stress  occurring  in  the  upper  surface  of  the  skin  for  the 
xmcracked,  patched  configuration  to  the  yield  stress  value,  corresponding  to  a  3.6  mm 
skin,  indicates  that  there  may  be  a  case  for  patching  the  complete  fleet  of  F-111  aircraft. 

Comparisons  made  between  simple  lap  joint  formulae  indicate  that  secondary  bending 
can  increase  the  adhesive  shear  and  peel  stresses.  However,  the  bending  correction 
factors  used  in  the  original  design  were  adequate.  The  differences  in  m.s.  between  the 
3D  F.E.  and  the  closed  form  solutions  are  not  significant. 

While  a  comparison  of  the  Von  Mises  adhesive  stresses  at  the  ends  of  the  patch  has 
shown  that  a  configuration  with  square  comers  may  have  advantages,  a  high  M.S. 
exists  at  this  location. 

The  Von  Mises  stresses  in  the  skin  are  lower  with  the  square  comer  configuration, 
however  this  is  not  significant. 

Not  only  has  3D  F.E.  confirmed  the  stractxual  integrity  of  the  RAAF  designed 
boron/ epoxy  patch  at  DUL,  but  also  the  use  of  those  equations  in  C5033  at  the  time  of 
the  design. 
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model  has  been  validated  using  strain  gauge  results  from  a  fuU  scale  test  wing.  The  results  of  the  F.E.  analysis  are 
shown  to  compare  favourably  with  closed  form  solutions  used  by  the  RAAF  (  RAAF  Engineering  Standard  C5033  ) 
in  the  original  design  of  the  repair.  Thus  the  present  work  provides  a  basis  for  confidence  in  the  design  procedures 
contained  in  RAAF  Engineering  Standard  C5033. 
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